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In  this  thesis,  we  propose  a  timing  model  for  fault  simulation  for  large  MOS  circuits 
which  effectively  balances  the  need  for  speed  and  accuracy.  Like  RSIM,  this  simulation 
model  represents  a  transistor  with  a  linear  resistor  and  uses  a  simple  RC  approximatior 
for  timing  estimates,  but  we  achieve  greater  accuracy  through  the  use  of  five  logic 
values  and  an  improved  timing  model.  Related  timing  simulators  have  been  developed 
by  other  researches.  In  addition,  our  simulator  can  handle  faults  such  as  stuck 
faults,  resistive  shorts,  and  threshold,  voltage  shifts.  This  is  possible  because 
our  timing  model  can  accurately  predict  the  behavior  of  the  circuit  under  these 
faults. 
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Due  to  its  rapid  and  accurate  evaluation  of  node  state  and  delays,  our  timing  model 
can  be  used  to  simulate  large  MOS  circuits.  In  addition,  it  can  be  easily  implemented 
as  a  concurrent  fault  simulator  which  w  ould  reduce  total  simulation  time  even  further.. 


The  remainder  of  this  thesis  is  organized  as  follows.  In  Chapter  2,  the  proposed 
timing  model  and  the  fault  model  will  be  described.  The  implementation  details 
of  FACT,  Fault  Simulation  with  Accurate  timing,  will  be  discussed  in  Chapter  3,  and 
in  Chapter  4,  results  on  various  circuits  will  be  shown  and  compared  with  those  of 
SPICE  and  RSIM. 
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CHAPI ER  1. 


INTRODUCTION 


With  the  increasing  need  for  manufacturers  to  maintain  quality  requirements  for  high  performance  and  high 
density  MOS  VLSI  integrated  circuits  and  the  multitude  of  physical  failures  and  defects  that  can  occur  in  such  cir¬ 
cuits,  delay  fault  testing  is  gaining  importance  [1,2,3. 4,5].  An  extremely  important  component  of  such  a  delay 
test  generation  environment  is  an  accurate  yet  fast  delay  fault  simulator.  The  goal  of  fault  simulation  is  to  gen¬ 
erate  the  responses  of  digital  circuits  under  both  fault  and  fault-free  conditions,  without  prohibitive  cost  Most 
simulators  available  today  either  sacrifice  accuracy  or  are  very  time-consuming.  Therefore,  choosing  the  right 
simulation  model  is  the  key  in  fault  simulation. 

MOS  circuit  models  that  were  previously  developed  for  true  value  simulation,  after  being  modified  to 
accommodate  fault  models,  have  been  used  by  fault  simulators,  such  as  FMOSSIM  [6].  These  fault  simulators  are 
incapable  of  detecting  timing  errors  and  even  some  logic  errors,  both  of  which  occur  in  actual  failures  in  the  MOS 
circuits  [7].  Others  have  shown  that  some  faults  only  affect  the  circuit  timing  and  the  delay  may  invalidate  test 
sets  [8,9].  Therefore,  fault  simulators  must  consider  the  circuit  delay  in  order  to  evaluate  test  sets  accurately. 
FAUST  [10]  is  one  such  fault  simulator  which  very  accurately  detects  timing  as  well  as  logical  errors.  It  was 
based  on  using  a  combination  of  table  lookup  and  numerical  integration  techniques  to  solve  differential  equaltions. 
With  today’s  large  MOS  circuits,  it  is  desirable  to  propose  faster  timing  models  without  sacrificing  the  accuracy 
too  much,  to  complete  the  simulation. 

In  this  thesis,  we  propose  a  timing  model  for  fault  simulation  for  large  MOS  circuits  which  effectively  bal¬ 
ances  the  need  for  speed  and  accuracy.  Like  RSIM  [II],  this  simulation  model  represents  a  transistor  with  a  linear 
resistor  and  uses  a  simple  RC  approximation  for  timing  estimates,  but  we  achieve  greater  accuracy  through  the  use 
of  five  logic  values  and  an  improv'd  timing  model.  Related  timing  simulators  have  been  developed  by  other 
researches  [12,13,14],  In  addition,  our  simulator  can  handle  faults  such  as  stuck  faults,  resistive  shorts,  and 
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threshold  voltage  shifts.  This  is  possible  because  our  timing  model  can  accurately  predict  the  behavior  of  the  cir¬ 
cuit  under  these  faults. 

Due  to  its  rapid  and  accurate  evaluation  of  node  state  and  delays,  our  timing  model  can  be  used  to  simulate 
large  MOS  circuits.  In  addition,  it  can  be  easily  implemented  as  a  concurrent  fault  simulator  which  would  reduce 
total  simulation  time  even  further. 

The  remainder  of  this  thesis  is  organized  as  follows.  In  Chapter  2,  the  proposed  timing  model  and  the  fault 
model  will  be  described.  The  implementation  details  of  FACT,  Fault  Simulation  with  ACcurate  Timing,  will  be 
discussed  in  Chapter  3,  and  in  Chapter  4,  results  on  various  circuits  will  be  shown  and  compared  with  those  of 
SPICE  and  RSIM. 
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CHAPTER  2. 

TIMING  MODEL  FOR  FAULT  SIMULATION 


A  new  model  for  timing  simulation  is  proposed  in  this  chapter.  The  first  section  describes  the  transistor 
model.  Using  this  model,  an  MOS  network  is  simulated  as  a  resistor  network  where  each  node’s  value  is  deter¬ 
mined  by  the  resistance  of  its  connections  to  various  inputs.  The  second  section  describes  the  node  model.  This 
is  followed  by  an  explanation  of  the  method  for  calculating  the  value  of  each  node.  The  fourth  section  discusses 
the  use  of  the  model  to  predict  the  propagation  of  new  input  values  through  a  network.  The  fifth  section  shows 
the  use  of  an  input  slope  and  an  output  load  capacitance  for  more  accurate  timing  predictions,  and  the  last  section 
describes  the  fault  model. 

2.1.  FACT  Model  of  Transistor 

In  our  delay  fault  simulator.  Fault  simulation  with  ACcurate  Timing  (FACT),  a  transistor  is  modeled  by  a 
linear  resistor  with  a  switch  in  series.  To  see  how  the  model  works,  let’s  consider  a  simple  inverter  in  Figure  2.1. 
We  can  think  of  the  effective  resistance  of  its  transistor  as 


The  figure  shows  the  actual  effective  resistance  of  both  pullup  and  pulldown  transistors  as  a  function  of  the 
inverter’s  output  voltage.  We  can  see  that  the  effective  resistances  of  the  uansislOio  change  as  uheir  terminal  vol¬ 
tages  vary.  However,  "average  channel  resistances"  might  be  enough  to  characterize  the  u-ansistor’s  behavior. 

Another  important  characteristic  of  the  transistor  is  that  it  behaves  like  a  switch.  With  certain  voltages 
applied  at  a  transistor’s  terminal  nodes,  it  makes  no  connection  between  its  source  and  drain  terminals  -  the 
transistor  is  "off."  As  the  terminal  voltages  change,  the  transistor  turns  "on,"  conducting  current  between  its 


source  and  drain  terminals. 
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Vdi:pulldo  wn 

Figure  2. 1  Effective  resistances  in  an  inverter 


In  Figure  2.2  are  the  three  basic  types  of  the  transistor  found  in  MOS  circuits.  The  only  difference  between 
the  n-channel  and  p-channel  transistors  is  the  gate  voltage  level  which  uutis  on  the  transistor.  The  depletion 
transistor  is  usually  connected  to  VDD  and  used  with  its  switch  always  on. 
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ON  when  gate  <  4V 

always  ON 

OFF  when  gate  <  IV 

OFF  when  gate  >  4V 

(a)  n-chaiuiel 
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switch 
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Figure  2.2.  Three  types  of  MOS  transistors 

From  the‘'v;  cjscrvations  on  the  MOS  transistor,  one  can  consinjcl  the  following  transistor  model  that  FACT 


uses: 


drain 


source 


Rtf 


source 


open  Vg„e  >  1 V 
closed  Vg.te  <  1 V 


(a)  n-channel 
transistor 


(b)  FACT  model 


Figure  2.3.  FACT  model  for  n-channel  transistor 


drain 


T 


source 


open  Vgjie  >  4V 
closed  Vg.ie  <  4V 


drain 

( I 

source 


(a)  p-channel 
transistor 
model 


(b)  depletion 
transistor 
model 


Figure  2.4.  FACT  models  for  p-channel  and  depletion  transistor 


Since  R^ff  is  the  only  device  information  FACT  uses  about  a  transistor  as  seen  in  the  above  figures,  choosing  the 
value  for  Retr  is  very  important  in  our  simplistic  Uansistor  model.  In  FACT,  the  effective  resistance  Res  is  deter 
mined  separately  for  each  transistor  and  depends  on 
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width,  length  Dimensions  of  the  active  transistor  area, 

type  Most  MOS  circuits  contain  more  than  one  type  of  transistor, 

context  Most  MOS  transistors  are  used  in  various  contexts, 

input  slope  Slope  with  which  the  input  changes  its  value. 

output  load  The  load  capacitance  at  the  output  that  needs  to  be  either  charged  or  discharged, 

gate  voltage  level  This  determines  how  "on"  the  transistor  is. 

Among  the  above  parameters,  it  is  the  gate  voltage  level  that  affects  the  value  of  R^fr  the  most.  Therefore,  Res  is 
first  determined  by  the  transistor’s  gate  voltage  level.  To  see  how  this  is  done,  consider  Figure  2.5  which  shows 
V-I  characteristics  of  an  MOS  transistor  with  various  voltages  applied  at  the  gate. 

In  the  figure,  as  the  gate  voltage  changes  from  one  to  five  volts,  the  V-1  curve  of  the  transistor  also  varies. 
Because  we  use  five  logic  levels  to  represent  a  node  state  (this  will  be  explained  in  the  following  section),  we  have 
a  V-I  curve  corresponding  to  each  of  five  logic  levels  applied  at  the  gate.  Modeling  a  tfansistor  by  a  linear  resistor 
is  equivalent  to  replacing  each  of  the  above  V-I  curves  with  a  straight  line,  with  its  slope  equal  to  the  conductance 
of  Reff 

The  values  for  R«ff,  the  inverse  of  the  slope  of  each  linearized  V-1  curve,  must  be  chosen  such  that  it  leads  to 
accurate  prerlictions  for  both  the  final  value  and  transition  time.  In  FACT,  this  is  done  in  a  pre-simulation  phase, 
in  which  a  series  of  experiments  are  performed  to  measure  the  resistance  of  the  transistor  in  various  circuit  con- 


Figure  2.5.  V-I  characteristics  of  an  MOS  transistor 


7 


texts.  Ideally,  the  experiments  should  be  performed  using  actual  circuits,  but  an  accurate  circuit  simulator  such  as 
SPICE  can  be  used  to  gather  the  needed  measurements.  At  this  point,  we  would  like  to  point  out  that  Re?,  actu¬ 
ally,  consists  of  three  resistances; 

R  niQc  when  calculating  the  final  voltage. 

R  dyniow  whcn  Calculating  the  transition  time  for  high-to-low  transitions. 

R  dynhigh  when  calculating  the  transition  time  for  low-to-high  transitions. 

Three  resistances  are  used  because  one  resistance  value  cannot  give  accurate  predictions  for  both  the  final  value 

and  transition  time. 

Static  resistances,  that  are  used  to  estimate  node  voltages,  are  comparatively  easy  to  choose.  When  all  the 
nodes  in  a  stage  are  connected  to  only  one  polarity  of  input,  the  values  chosen  for  static  resistances  do  iiOt  affect 
the  voltage  computadon.  When  a  circuit  makes  a  connection  to  inputs  of  different  polarities,  the  transistor  resis¬ 
tances  must  be  chosen  to  predict  the  output  voltage.  Since  only  the  ratio  of  the  puUup  and  pulldown  devices  is 
constrained,  there  is  considerable  freedom  in  choosing  the  actual  resistance  values. 

Choosing  the  appropriate  dynamic  resistance  is,  however,  not  so  simple.  Dynamic  resistances  are  used  to 
predict  the  transition  time.  Therefore,  each  of  the  experiments  consists  of  measuring  the  length  of  time  required  for 
the  output  to  rise  or  fall  from  its  starting  voltage  to  the  switching  threshold.  If  the  load  capacitance  is  known, 
Rdyn’jw  and  Rdynhigh  c^n  be  calculated,  essentially  inverting  the  computation  performed  by  FACT.  As  an  example, 
consider  the  experiment  with  an  n-channel  pulldown  in  Figure  2.6.  Initially,  the  gate  is  at  zero  volt  and  the  capa¬ 
citance  at  the  output  node  is  fully  charged  to  5  volts.  By  applying  various  voltages  to  the  gate,  we  turn  the 
transistor  on  and  discharge  the  output  capacitance.  This  is  followed  by  observing  the  output  waveform  and 
approximating  it  with  an  exponential  curve.  We,  then,  measure  the  delay,  t ,  the  time  taken  for  the  output  voltage 
to  ''  op  to  36.7%  of  its  initial  value,  the  transition  threshold,  and  obtain  Rdyniow  which  is  equal  to  The  values 
currently  used  in  FACT  are  given  in  Figure  2.7. 

Once  obtained  in  the  pre-simulation  phase,  Refr  are  used  throughout  the  actual  simulation  phase.  However, 
as  seen  in  the  previous  example,  the  method  for  obtaining  R«£f  does  not  take  the  input  slope,  the  load  size,  and  the 
transistor  size  into  account  Therefore,  the  effective  resistance  must  go  through  a  series  of  calibrations  taking 
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Figure  2.6.  An  experiment  with  an  n-channel  pulldown  to  measure  R^r 
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n-mos 

0->2V 
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315 

0-»3V 

130 

80 

0-»4V 

65 

35 

0^5V 

35 

21 

p-mos 

5^3V 

315 

510 

5  ->2V 

80 

130 

5  ^  IV 

35 

65 

5  -♦OV 

21 

35 

Figure  2.7.  Rjynimic  used  by  FACT 


these  parameters  into  account  before  it  is  ready  to  be  used.  This  calibration  phase  will  be  discussed  in  detail  in 
Section  2.4. 


2.2.  FACT’S  Node  Model 


To  simulate  a  digital  circuit,  at  least  two  logic  states,  logic-one  and  logic-zero,  are  necessary  to  represent  the 
node  voltage.  Also,  a  third  logic  state,  referred  to  as  the  "X”  stale,  is  often  used  to  represent  the  unknown  or  the 
intermediate  voltage  levels.  Representing  the  node  by  this  binary  or  ternary  logic  values  may  be  sufficient  if  the 


circuit  is  designed  so  conservatively  that  intermediate  voltage  levels  rarely  occur.  Unfortunately,  in  today’s  digital 
MOS  circuits,  intermediate  voltage  levels  fifequently  occur.  For  example,  let  us  consider  the  circuit  in  Figure  2.8. 
In  the  figure,  the  gate  voltage  at  node  A  turns  the  transistor  T1  on  to  charge  up  the  capacitance  at  node  B.  How¬ 
ever,  the  voltage  at  node  B  rises  only  up  to  4  volts,  instead  of  the  desired  voltage,  5  volts,  due  to  the  threshold 
voltage  drop  between  the  gate  and  the  drain  node  of  transistor  Tl.  This  weak  logic-one  value  is,  next,  used  to  turn 
the  next  transistor  T2  on.  Due  to  the  threshold  voltage  drop  again,  the  voltage  at  node  C  is  pulled  only  up  to  3 
volts,  an  intermediate  voltage. 

Intermediate  voltages  can  also  arise  by  a  charge  sharing  which  happens  when  two  source  free  transistor  net¬ 
works  become  one  by  a  source  and  a  drain  connection  of  a  transistor.  As  soon  as  the  connection  is  made,  nodes  in 
the  network  are  allowed  to  share  charge  and  all  will  reach  the  same  intermediate  voltage  level. 

Intermediate  voltage  levels  also  occur  when  a  fault  is  injected  into  the  circuit  Bryant  of  FMOSSIM  [6] 
showed  that  many  fault  models  spread  intermediate  voltages  throughout  the  circuit.  For  an  accurate  fault  simula¬ 
tion,  therefore,  intermediate  voltages  should  be  modeled  accurately. 

Modeling  intermediate  node  voltages  is  also  important  because  it  directly  affects  the  operation  of  the  transis- 
tw,  it  is  responsible  for  the  different  ''on"  sf.tes  of  a  transistor.  For  example,  if  the  gate  of  an  NMOS  transistor  is 
driven  by  a  full  5  volts,  the  transistor  is  fully  "on"  and  behaves  like  a  closed  switch.  As  the  gate  voltage 
decreases,  however,  the  transistor  starts  acting  like  a  resistor,  with  its  resistance  increasing  as  the  gale  voltage 

OV  5V 

X. 

4V 

0V-^3V 

c 

Figure  2.8.  A  circuit  showing  an  occurrence  of  intermediate  voltages 
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decreases.  Finally,  it  becomes  an  open  switch  when  the  gate  voltage  drq)s  to  and  below  one  volt. 

In  (xder  to  build  an  accurate  fault  simulator,  we  chose  to  quantize  the  node  into  five  logic  levels  (See  Figure 
2.9).  Now,  we  can  not  only  describe  the  transistor’s  behavior  more  accurately  but  also  model  the  intermediate 
voltages  firequently  occurring  in  fault  simulation.  A  detailed  study  on  the  effect  of  physical  failure  of  MOS  cir- 
cuits[lS]  has  also  shown  that  five  logic  levels  are  necessary  to  represent  the  intermediate  voltages  caused  by  a 
wide  variety  of  fault  models. 

2J.  FACT’S  Final  Value  Computation 

The  final  value  of  a  node  is  obtained  as  follows.  Using  the  transistor  model  described  in  the  previous  sec¬ 
tion,  the  original  network  is  transformed  into  a  network  of  resistors  and  capacitors  as  seen  in  Figure  2.10  and  Fig¬ 
ure  2.1 1.  When  a  node  evaluation  takes  place,  the  evaluation  is  performed  not  on  the  node  in  isolation,  but  on  an 
entire  set  of  nodes  that  are  cunendy  connected  via  source-to-drain  connections  of  non-off  transistors.  We  call 
such  a  group  of  nodes  a  stage  (See  Figure  2.12.).  Stages  do  not  extend  through  input,  power  or  ground  nodes. 

In  our  node  value  computation,  we  evaluate  the  stages  as  trees.  If  the  stage  has  any  cyclical  connections,  it 
is  nevertheless  analyzed  as  a  tree  by  arbitrarily  ignoring  any  transistors  which  complete  a  cycle  (See  Figure  2.13.). 
This  simplification  can  result  in  serious  errors  in  estimation;  fortunately,  stages  containing  cycles  are  rare. 


voltage  range 

logic  value 

0  —  IV 

0 

1  —  2V 

0* 

2  — 3V 

I 

3  — 4V 

1* 

4  — 5V 

1 

Figure  2.9.  Multiple  logic  values  used  by  FACT 


Figure  2.10.  Typical  MOS  Signal  Distribution  Network 


Figure  2.11.  Linear  circuit  model  for  the  network  in  Figure  2.10 


Figure  2.12.  Simple  circuit  that  has  three  stages 
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Vdd 


Figure  2.13.  Tree  approximation 

Now,  the  node  value  computation  is  equivalent  to  solving  a  system  of  linear  equations.  Quite  often,  the 
stage  has  only  one  driving  source,  in  which  case  the  solution  is  trivial.  But,  multiple  driving  sources  may  also 
occur  as  seen  in  Figure  2.14.  In  such  a  case,  FACT  solves  the  linear  system  exactly,  using  a  tree-based  algorithm. 

In  the  tree-based  algorithm,  given  a  particular  node,  the  efiFect  of  inputs  on  the  node  is  found  by  finding  the 
exact  Thevenin  equivalent  for  the  stage,  with  the  given  node  as  the  output  (See  Figure  2.15.). 


Figure  2.14.  A  multiple-source  tree 
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Figure  2.15.  Equivalent  circuit  for  a  network  node 

Values  for  all  the  nodes  in  a  stage  are  obtained  in  0(n)  by  passing  the  tree  twice.  For  example,  consider  the  RC 
tree  with  multiple  sources  in  Figure  2.16(a).  Since  we  are  only  interested  in  steady-state  node  values,  we  can 
ignore  capacitances  from  the  network.  The  algorithm,  first,  identifies  a  ‘root’  which  is  a  node  directly  connected 
to  a  source  through  a  resistor.  In  the  first  pass  of  the  tree,  Vuwv  and  Rthev  at  the  root,  looking  firom  the  source  side, 
ate  obtained  by  using  two  rules  shown  in  Figure  2. 17.  By  applying  these  rules,  we  can  reduce  the  original  networic 
into  the  equivalent  network  shown  in  Figure  2.16(b).  The  exact  voltage  at  the  root  node  is,  then,  found  by  apply¬ 
ing  Rule  2.  The  rest  of  the  node  values  arc  obtained  in  the  second  pass  by  walking  back  to  the  tree  network. 

In  the  second  pass,  Vihev  and  Rihcv,  looking  toward  the  root  node,  are  found  for  each  node  in  the  network  as 
seen  in  Figure  2.16(c).  At  each  node,  these  Thevenin  equivalent  values  are  combined  with  those  obtained  in  the 
first  pass  to  obtain  the  node  value  by  using  Rule  2.  The  tree  walk  in  the  second  pass  continues  to  proceed  until  a 
source  node  is  reached.  A  similar  approach  was  taken  by  VTIsim  [16]. 

2.4.  FACT’S  Transition  Time  Computation 

The  aim  of  FACT  is  to  generate  the  time  response  of  digital  circuit  designs  with  less  detail  and  at  higher 
speed  than  circuit  level  simulation,  but  without  losing  accuracy  as  measured  in  terms  of  voltage  levels  and  timing. 
This  is  achieved  by  modeling  the  nonlinear  transistor  network  as  linear  resistor  networks  as  previously  described, 
and  apply  an  RC  network  technique  to  find  the  delay. 
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R2 


(a)  An  RC-tree  network 


Thevenin 

equivalent 


(b)  Thevenin  equivalent  in  the  first  pass 


Thevenin 

equivalent 


(c)  Thevenin  equivalent  in  the  second  pass 


Figure  2.16.  Tree-based  algorithm  for  final  value  compuiauon 
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(a)  Rule  1. 


RrR2 


(b)  Rule  2. 

Figure  2.17.  Rules  to  obtain  the  Thevenin  equivalent  circuit 

FACT’S  method  for  delay  computation  is  based  on  the  Penfield-Rubinstein’s  models  for  delay  estimate  [17]. 
Penfield  and  Rubinstein  derived  a  tight  upper  ai  d  lower  bound  for  a  delay  waveform  in  an  RC  tree  and  proposed  a 
delay  estimate  that  always  lies  within  the  bounds.  To  see  how  the  model  works,  consider  the  RC  tree  in  Figure 
2.18. 

An  RC  tree  is  a  resistor  tree  where  the  root  of  the  tree  is  the  input  and  each  node  in  the  tree  has  a  capacitor 
to  ground.  Since  the  only  dc  path  to  a  node  is  from  the  input,  given  enough  time  all  the  nodes  in  the  tree  will  set¬ 
tle  to  the  same  voltage:  the  voltage  at  the  input  node.  When  a  step  voltage  is  applied  to  the  input,  all  the  nodes  in 
the  tree  lag  the  input,  eventually  settling  to  the  new  value.  Because  the  output  waveforms  change  gradually,  the 
most  complete  method  of  specifying  the  delay  for  an  input  voltage  step  is  to  find  the  output  waveforms  versus 
time.  The  delay  estimate  proposed  by  Penfield  and  Rubinstein  trys  to  approximates  this  output  waveform  with  an 


exponential. 


16 


Consider  any  output  node  e,  and  any  capacitor  at  node  k  with  capacitance  Ck.  The  resistance  Rke  is  defined 
as  the  resistance  of  the  portion  of  the  unique  path  between  the  input  and  e,  that  is  common  with  the  path  between 
the  input  and  node  k.  The  approximate  output  waveform  at  node  e  is  defined  as 

Ve(t)  =  e'^;  tDt  =  ^  RckCk 

This  simple  delay  estimate  takes  the  parasitic  capacitances  at  all  other  nodes  in  the  stage  into  account  This 
method  is,  therefore,  likely  to  produce  more  accurate  results  than  those  of  table-driven  methods  [18,19,20],  Its 
efficiency  in  computation  also  makes  it  suitable  for  a  fast  timing  simulator  such  as  FACT.  One  disadvantage  of 
this  delay  model  is  that  it  is  only  applicable  to  an  RC  tree:  it  cannot  be  used  for  circuits  with  multiple  driving  vol¬ 
tages.  Because  many  sub-circuits  have  multiple  same-polarity  source  connections,  crude  approximations  are  una¬ 
voidable  if  this  model  is  used  for  such  circuits.  For  an  accurate  simulation,  therefore,  it  is  necessary  to  extend  this 
model  to  handle  non-RC  trees. 

In  FACT,  this  is  achieved  by  a  computationally  efficient  algorithm  we  developed  based  on  Wyatts’s  work  on 
delay  estimates  in  RC  meshes  [21],  Wyatt  has  shown,  with  redefinition  of  Rke,  that  the  exponential  estimate  for  an 
RC  tree  can  also  be  applied  to  an  RC  mesh.  In  the  RC  mesh  in  Figure  2.19,  the  voltage  at  each  node  can  be 
represented  by  the  following  matrix,  with  a  new  definition  for  rke 
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R2  R3 


vi-e 

rii  ri2  ri3  ri4 

ii 

vr-e 

til  r22  r23  r24 

i2 

V3-e 

r3i  r32  r33  r34 

is 

y4-e 

r4i  r42  r43  r44. 

-U. 

fke  =  where  ij  =  0  for  all  j>0,  j/k 

Unlike  the  Rke  of  RC  trees,  which  can  be  read  off  by  observing  the  circuit,  calculation  of  rke  is  quite  com¬ 
plex.  For  example,  iv*  can  be  obtained  by  inverting  an  (n-'l)x(n-l)  conductance  matrix,  where  n  is  the  number  of 
non-ground  nodes  in  circuit  This  matrix  will  have  positive  diagonal  elements  and  non-positive  off-diagonal  ele¬ 
ments.  Each  node  except  for  ground  corresponds  to  one  row  and  one  column  of  the  conductance  mauix.  Entries 
in  the  matrix  are  determined  as  follows.  The  diagonal  entry  of  row  j  is  equal  to  the  sum  of  the  conductances 
incident  on  node  k.  The  off-diagonal  entries  are  zero  unless  node  i  corresponding  to  that  column  is  a  neighbor  of 
node  k.  In  this  case  that  entry  is  set  equal  to  the  negative  of  the  conductance  between  i  and  k.  It  is  easily  shown 
that  this  is  equivalent  to  writing  KCL  at  node  k.  The  resulting  conductance  matrix  is  augmented  with  an 
(n-l)x(n-l)  identity  matrix.  Gaussian  elimination  or  the  Grout  algorithm  [22]  is  used  to  reduce  the  conductance 
matrix  to  an  identity  matrix.  At  this  point,  the  augmented  matrix  holds  the  inverted  conductance  matrix.  Thus, 
[0:1]  -»  [I:  G"'].  The  G"'  matrix  contains  the  resistance  values  defined  by  Wyatt.  Since  the  inversion  of  a  matrix 


requires  an  0(  n^ )  algorithm,  this  method  is  computationally  too  expensive. 
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For  an  RC  tree  with  a  multiple  same-polarity  source  connection,  it  is  possible  to  obtain  efficiently 
without  having  to  invert  a  matrix.  This  is  done  by  using  a  tree-based  algorithm  we  developed.  To  see  how  the 
algorithm  works,  consider  the  RC  tree  driven  by  two  sources  in  Figure  2.20.  As  seen  in  the  figure,  we  partition 
the  tree  with  multiple  sources  into  sub-trees  such  that  the  root  of  each  sub-tree  is  directly  connected  to  each  source 
through  a  resistor.  By  representing  an  RC  tree  as  an  RC  forest  in  this  way,  we  can  easily  find  the  necessary  resis¬ 
tances,  rke.  First,  consider  the  case  in  which  both  nodes  e  and  k  belong  to  the  same  sub-tree.  For  the  purpose  of 
illustration,  let  k  and  e  be  any  arbitrary  nodes  in  Tree  1  of  Figure  2.18;  k  =  G  and  e  =  H.  By  definition, 

roH  =  -^2.  with  all  the  currents  across  the  capacitor  at  each  node  set  to  zero  with  the  exception  of  in  .  Therefore. 
tH 

in  originating  from  the  capacitor  at  node  H  can  flow  to  ground  only  through  the  resistive  paths, 
R9->R5-»R3-+R2  and  R9->R5— >R1.  All  other  resistors  can  be  ignored  from  the  picture  because  no  current  flows 
through  them.  The  resulting  circuit  is  shown  in  Figure  2.21. 


Since  the  current  through  R8  is  also  zero,  Vo  =  Vo  by  the  KVL  and  b  =  b  by  the  KCL.  Rgh  is,  therefore, 

equal  to  which  is  the  equivalent  resistance  between  node  D  and  ground.  The  problem  of  finding  Rke  reduces 
b 


to  finding  this  equivalent  resistance  at  node  D.  Rather  than  trying  to  find  this  resistance  directly,  we  divide  this 


Figure  2.20.  Multiple  source  tree 
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equivalent  resistance  into  two  parts;  Riourc<i  the  resistance  between  the  root  node  and  the  source  nodc>  and  Riiee* 
the  resistance  of  the  unique  path  between  the  root  and  e,  that  is  common  w-th  the  path  between  the  root  and  node 
k.  Rke  is  simply  the  sum  of  R  sub  source  and  R  sub  tree. 

The  computation  method  for  Rtm  is  identical  to  the  one  used  in  calculating  Rke  of  a  one-source  RC  tree. 
Therefore,  we  can  treat  the  n-source  RC  tree  as  n  separate  one-source  RC  trees  and  use  the  same  algorithm  If 
nodes  e  and  k  are  not  in  the  same  sub-tree,  the  path  from  node  e  to  the  root  of  its  sub-tree  is  exclusive  from  the 
path  from  node  k  to  die  root.  Therefore,  Rk*  is,  simply,  equal  to  Rwun*  in  this  case. 

After  obtaining  Rut*,  we  compute  Riouce  and  add  it  to  Rue*.  Computing  Ruw  requires  an  0(n)  algorithm, 
where  n  is  the  number  of  nodes  in  the  stage,  and  RKwree  is  obtained  by  using  an  0(s^)  algorithm,  where  s  is  the 
number  of  sources.  Since  s  is  usually  limited  to  five  or  less  per  stage,  the  computation  is  dominated  by  the  first 
one.  Therefore,  the  complexity  of  the  above  algorithm  is  0(n). 

The  extension  of  the  delay  model  to  handle  an  RC  tree  with  multiple  source  enables  a  simulator  to  cover  a 
wide  variety  of  circuits  for  an  accurate  simulation.  For  instance,  delay  at  the  output  of  a  NOR  gate  can  be 
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acx:urately  computed  even  when  the  two  inputs  change  simultaneously.  The  effect  of  having  mu'.aple  paths  to  Vdd 
for  a  faster  pull-up  can  also  be  seen  clearly  in  our  timing  model. 

With  a  slight  modification,  our  timing  model  can  also  handle  dynamic  CMOS  logic  structures  accurately. 
For  example,  consider  the  basic  dynamic  CMOS  gate  in  Figure  2.22(a).  It  consists  of  an  n-transistor  logic  su^c- 
ture  whose  output  node  is  precharged  to  Vdd  by  a  p-transistor  and  conditionally  discharged  by  an  n-transistor  con¬ 
nected  to  Vj,.  4)  is  a  clock  signal.  The  precharge  phase  occurs  when  <|)  =  0.  The  path  to  V„  supply  is  closed  via 
the  n-transistor  when  0=1.  And,  if  the  inputs  to  the  n  logic  block  are  set  such  that  there  exists  a  path  from  the 
output  to  this  evaluation  transistor,  the  output  is  pulled  down  to  zero.  If  the  n  logic  block  docs  not  include  a 
cycle,  the  delay  at  the  output  can  be  accurately  modeled. 

In  its  evaluation  phase,  a  dynamic  logic  circuit  can  be  modeled  by  the  RC  circuit  in  Figure  2.22(b).  To  find 
the  delay  at  the  output,  we  need  to  find  Re,outpui  for  e  =  all  nodes  in  the  stage.  Let  c  be  an  arbitrary  node  in  the 
figure.  To  find  Re.ou(pui,  we  inject  current  lomput  at  the  output  and  measure  the  voltage  at  node  e.  We-  The  desired 
y 

ratio,  - — is  obtained  as  follows: 

Vdd  _^output 

LOG  id 
BLOCK,^ 

:  Vhlock 

V 

A 

i-f- 

|v„., 

Figure  2.22(a)  A  dynamic  CMOS  gate  (b)  An  equivalent  resistive  network 


21 


Ve  —  Vblock  Vcv*l 


SinC6  Vcval  —  loutpul’Rcvtl j 


Vc  —  Vwock  Iouipul*R«val 


Dividing  each  side  by  Ioutpui> 


=  .^i£i  + 

^output  ^tput 

If  Rev»i  =  0,  the  circuit  becomes  an  RC  tree  with  two  same-polarity  sources  and  Rc.ouiput'  for  this  RC  tree  equals 
■  Therefore,  R«,omput  for  the  dynamic  circuit  is  obtained  by  adding  Rev^  to  Re,ouiput'  of  this  RC  tree.  Re,outpui' 

'output 

is  obtained  by  using  the  same  algorithm  we  developed  for  multiple  source  RC  trees,  and  is  added  to  Rev.i; 

R«,output  “  Rc.output  ReviJ 

Therefore,  the  algorithm  for  multiple-source  RC  trees  can  also  be  used  for  such  dynamic  circuits. 


2S.  Slope  and  Load  Effect 

A  significant  source  of  error  in  the  RC  model  comes  from  its  inability  to  deal  with  waveform  shape.  In 
practice,  the  effective  resistance  of  a  transistor  depends  on  the  waveform  on  its  gate.  If  the  trigger  transistor  turns 
on  instantaneously,  then  its  full  driving  power  is  used  to  drain  the  output  capacitance  and  the  transistor  has  a  rela¬ 
tively  low  effective  resistance.  If  the  trigger  transistor  turns  on  slowly,  then  it  may  do  much  or  all  of  its  work 
while  only  partially  tumed-on.  In  this  case  its  effective  resistance  will  be  higher. 

If  all  waveforms  in  a  circuit  have  the  same  shape,  then  the  effective  resistances  of  transistors  can  be  charac¬ 
terized  using  that  waveform  and  the  RC  model  will  produce  accurate  results.  Unfortunately,  this  is  not  the  case  in 
actual  VLSI  circuits.  Although  almost  all  waveforms  have  an  exponential  shape,  they  vary  by  more  than  three 
orders  of  magnitude  in  their  slopes.  As  a  result,  the  effective  resistance  of  the  transistors  varies  by  more  than  a 
factor  of  ten.  Therefore,  it  is  highly  desired  to  include  the  effect  of  the  input  slope  in  the  calculation  of  the 


effective  resistance. 
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Unfortunately,  the  effective  resistance  of  a  transistor  depends  not  only  on  the  slope  of  its  gate  voltage,  but 
also  on  the  load  being  driven  by  the  stage  and  on  the  sizes  of  the  transistors  in  the  stage.  If  a  stage  is  driving  a 
large  load,  or  has  very  small  transistors,  then  only  very  slowly  rising  input  slopes  will  affect  the  stage’s  delay.  If  a 
stage  is  driving  a  small  load  or  has  very  large  transistors,  its  delay  will  be  more  sensitive  to  the  slope  of  its  input. 

As  shown  previously,  effective  resistances  are  obtained  with  the  assumption  that  the  input  is  driven  by  a  step 
function.  Therefore,  if  the  input  is  actually  a  step  voltage,  this  value  will  be  used  as  its  effective  resistance.  For 
inputs  with  slopes.  Re?  is  modified  as  follows: 

Here,  k  is  a  fudge  factor  that  is  determined  in  the  pre-simulation  phase  by  using  SPICE.  W  and  L  are  the  size  of 
the  transistor  in  X.  The  load  is  approximated  by  the  sum  of  capacitances  driven  by  the  transistor. 

The  accuracy  of  this  model  depends  largely  on  the  accuracy  with  which  slope  can  be  calculated.  The 
current  FACT  implementation  approximates  the  slope  by  using  the  delay  of  the  input.  For  any  given  voltage,  the 
slope  of  an  exponential  waveform  at  that  voltage  is  proportional  to  the  delay  time  to  reach  that  voltage.  This 
means  that  the  slope  at  a  node  is  proportional  to  its  delay,  which  is  the  delay  computed  by  the  RC  timing  model. 
A  similar  approach  was  also  taken  by  the  Crystal  timing  analyzer  [23]. 

2.6.  Fault  Model  of  FACT 

The  faults  that  can  be  modeled  in  FACT  are:  (1)  Node  stuck-at-zero  (2)  Node  stuck-at-one  (3)  Short 
between  nodes  (4)  Line  open  (5)  Transistor  stuck-open  (6)  Transistor  stuck-closed  (7)  Gate-to-drain,  gate-to- 
source,  source-to-drain  shorts  of  transistors.  (8)  Threshold  voltage  shifts  of  transistors  giving  rise  to  changes  in 
characteristics.  In  FACT,  each  of  the  shorts  and  opens  described  can  be  "resistive"  with  various  ranges  of  resis¬ 
tance  values.  This  is  possible  because  these  faults  are  represented  by  adding  extra  fault  transistors  such  as  FMOS- 
SIM  [6],  and  any  desired  resistance  of  the  fault  can  be  obtained  by  adjusting  the  gate  voltage  and  the  width  of  the 
fault  transistor. 


Figure  2.23  shows  various  faults  and  the  corresponding  model.  Those  gate  nodes  labeled  f  are  normally  0, 
but  are  set  to  0*.  I,  I*,  and  1  to  activate  the  fault  with  various  resistances;  gate  nodes  labeled  f  are  normally  1,  but 
are  set  to  0,  0*,  I,  and  I*  to  activate  the  fault.  Those  gate  nodes  labeled  s  are  normally  0,  but  are  set  to  1  to  short 
the  source  and  the  drain  of  the  transistor;  gate  nodes  labeled  s  are  normally  1,  but  are  set  to  0  to  open  the  source 
and  the  drain  nodes.  The  effective  resistances  of  these  transistors  are  set  to  0  ohm  to  have  the  effect  of  a  short 
between  the  source  and  the  drain  nodes. 

A  stuck-at-zero  or  stuck-at-one  node  fault  can  be  modeled  by  inserting  a  transistor  with  an  infinite  width  to 
short  the  node  to  Gnd  or  Vdd,  respectively.  A  stuck-ciosed  transistor  fault  is  injected  by  shorting  the  transistO'’’s 
source  and  drain  nodes.  The  short  is  achieved  by  a  fault  transistor  in  parallel  with  the  fault-free  transistor.  To 
eliminate  the  effect  of  the  fault-free  transistor  when  simulating  the  fault,  we  insert  an  S-transistor,  which  acts  like 
a  switch,  at  the  source  node  of  the  fault-free  transistor.  Similarly,  a  stuck-open  transistor  fault  is  modeled  by  a 
fault  transistor  in  series  and  an  S-transistor  in  parallel  as  seen  in  Figure  2.23.  A  gate-to-source  short  of  a  transistor 
is  injected  by  connecting  a  fault  uansistor  to  the  gate  and  the  source  terminals  of  the  transistor.  Gate-to-drain  and 
source-to-drain  shorts  are  modeled  similarly.  Threshold  voltage  shifts  of  transistors  are  injected  by  shifting  the 
values  for  R«e  in  discrete  steps. 
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Figure  2.23.  Fault  models  used  by  FACT 
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CHAPTER  3. 

IMPLEMENTATION 


FACT  uses  the  event-driven  method,  which  essentially  propagates  the  effects  of  signal  changes  through  a 
circuit  The  main  operations  to  be  performed  in  an  event-driven  simulation  algorithm  are  scheduling  of  events, 
processing  of  events,  evaluation,  and  event  list  management. 

3.1.  Scheduling 

Events  are  generated  as  a  result  of  primary  input  changes  and  the  evaluation  of  stages.  All  events  for  a  par¬ 
ticular  time  are  kept  in  a  list,  called  the  event  list,  for  that  time.  Insertion  of  events  in  the  event  list  for  the 
appropriate  time  is  termed  scheduling. 

Only  signal  changes  are  inserted  into  the  event  list.  This  requires  a  comparison  of  the  computed  value  of 
each  node  with  its  present  value  if  a  change  has  not  already  been  scheduled.  If  a  change  has  been  scheduled,  the 
newly  computed  value  is  compared  with  the  last  scheduled  value.  Suppression  of  spikes  of  short  duration  is  per¬ 
formed  at  this  time. 

For  detecting  a  spike,  the  newly  computed  value  of  a  signal  is  compared  with  its  last  scheduled  value.  A 
spike  is  present  if  the  values  are  different  and  the  time  between  the  last  scheduled  change  and  the  new  change  to 
be  scheduled  is  smaller  than  a  prescribed  value.  The  spike  is  suppressed  by  deleting  the  last  scheduled  event  and 
ignoring  the  new  change. 

In  FACT,  it  is  sometimes  necessary  to  cancel  certain  scheduled  events  in  order  to  obtain  more  accurate  uan- 
sition  time.  For  example,  consider  a  CMOS  NOR  gate  in  Figure  3.1,  where  input  IN2  rises  immediately  after 
input  INI  does.  When  input  INI  changes  its  state,  the  output  gradually  goes  to  zero,  and  a  new  event  is  scheduled 
at  the  transition  time,  x.  But  before  x  is  reached,  input  IN2  changes  its  state  and  turns  on  an  n-channel  transistor, 
adding  another  pulldown  to  the  output.  A  new  uxmsition  time  is  computed  using  the  two  pull-downs  and  is 
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Figure  3.1.  Effect  of  multiple  changing  inputs 


scheduled  before  t.  The  previously  scheduled  event  at  t  should,  therefore,  be  canceled  and  replaced  by  the  newly 
computed  one. 

3.2.  Processing  Scheduled  Events  and  Stage  Evaluation 

The  processing  of  scheduled  events  consists  of  updating  node  values  and  determining  the  stages  to  be 
evaluated  as  a  consequence  of  the  node  value  change.  For  every  event  in  the  event  list  at  the  current  simulated 
time,  the  node  value  is  updated  to  its  new  value.  If  the  updated  node  is  connected  to  gates  of  other  transistors,  a 
new  stage  is  formed  for  such  transistors  and  is  put  on  a  list,  called  the  evaluation  list,  if  it  is  not  already  there. 

After  all  the  updates  and  fan-out  processing  for  a  particular  simulated  time  have  been  performed,  the  stages 
on  the  evaluation  list  are  evaluated.  The  stage  evaluation  consists  of  computing  new  node  values  and  their  transi¬ 
tion  times.  Events  generated  as  a  result  of  the  stage  evaluation  are  inserted  into  the  event  list 

3J.  Event  List  Management 

FACT  uses  the  time  wheel  to  manage  the  event  list  (See  Figure  3.2.).  The  event  lists  on  the  time  wheel  are 
separated  by  fixed  increments  of  1  nanosecond  of  simulated  time,  allowing  faster  access  to  events  for  a  desired 
time.  The  event  wheel  is  managed  as  a  circular  buffer  in  which  the  N  array  elements  hold  events  for  the  next  N 
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Figure  3.2.  The  time  wheel  is  implemented  with  an  event  array 


time  quanta.  An  array  index  indicates  the  array  element  which  corresponds  to  the  current  simulated  time.  If  a 
new  event  is  scheduled  fw  a  time  M  quanta  in  the  future,  where  M  <  N,  the  event  is  added  to  the  end  of  the  event 
list  sorted  in  the  array  element  (index  +  M)  mod  N.  If  M  >  N,  the  event  is  inserted  into  the  overflow  list  accord¬ 
ing  to  its  scheduled  time. 

To  find  the  next  event  to  process,  the  event  array  is  searched  starting  at  the  current  index,  until  an  event  is 
found.  Each  increment  of  the  index  corresponds  to  advancing  the  simulated  time  by  one  nanosecond.  If  the  array 
is  empty,  the  simulated  time  is  advanced  to  equal  the  scheduled  time  of  the  first  event  on  the  overflow  list;  this 
event  then  becomes  the  next  one  to  be  processed.  When  an  event  is  located  for  processing,  the  overflow  list  is 
examined  to  find  events  whose  scheduled  times  are  less  than  N  time  quanta  away  from  the  new  simulated  time. 
Such  events  are  moved  from  the  overflow  list  to  the  appropriate  list  in  the  time  wheel. 
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CHAPTER  4. 

SIMULATION  RESULTS 


The  timing  simulation  model  described  in  Chapters  2  and  3  has  been  implemented.  The  fault  simulator  with 
an  accurate  timing  model,  FACT,  is  written  in  C  and  runs  on  the  SUN  workstation  under  the  Unix  operating  sys¬ 
tem.  We  evaluate  the  performance  of  FACT  based  on  its  computational  speed  and  the  accuracy  of  its  timing  esti¬ 
mates,  first,  for  some  circuits  without  fault  injection,  and  then  for  some  circuits  with  faults. 

In  Figure  4  i,  a  CMOS  NOR  gate  is  shown.  Here,  the  effect  of  the  second  input  change,  IN2,  is  clearly 
seen  by  the  change  in  the  output  slope  (Figure  4.2.).  This  is  possible  because  the  voltage  transitions  are  modeled 
by  ramps  in  FACT.  RSIMfll]  solves  this  event  scheduling  iroblem  by  ignoring  the  first  input  change,  resulting 
in  a  serious  timing  error. 

The  next  example  is  a  Manchester  carry  chain  (Figure  4.3.).  Here,  we  can  see  how  our  improved  delay 
algorithm  adds  accuracy  to  the  timing  estimate.  In  this  circuit,  the  length  of  delay  at  the  output  of  the  carry  chain 
depends  upon  the  number  and  the  positions  of  pull-down  paths.  This  variation  in  delay  is  more  accurately 
predicted  by  the  multiple-source  delay  algorithm  of  FACT  than  that  of  RSIM  (Figure  4.4). 

One  large  benefit  of  using  five  logic  values  over  0,  1  and  X  values  can  be  seen  when  simulating  circuits 
with  a  feedback  such  as  EXOR  gates  in  Figure  4.5.  Each  EXOR  gate  in  the  figure  contains  a  pass-gate  and  a 
feedback  line  which  controls  the  pass-gate.  For  certain  input  transitions,  the  logic  transition  at  the  output  has  an 
initial  slow  rise  time  due  to  a  slow  feedback  response.  RSIM  and  other  switch  level  simulators  fail  to  predict  the 
correct  value  at  the  output  because  the  small  voltage  change  at  the  output  is  not  recognized  as  an  event.  In  FACT, 
however,  the  correct  value  is  predicted  because  the  voltage  change  as  little  as  0.2V  is  recognized  as  an  event  (Fig¬ 
ure  4.6.).  This  sensitivity  added  to  the  event  scheduling  enables  the  simulator  to  predict  the  corrcct  output.  Table 
4.1  'hows  a  comparison  of  the  execution  times  of  RSIM,  FACT,  and  SPICE  for  various  circuits. 


29 


Figure  4.7  shows  a  chain  of  CMOS  inverters  with  a  gate-to-source  resistive  short  faulL  Figure  4.8(a)  shows 
the  output  waveform  at  node  D  when  R  =  200  ohms.  This  fault  causes  a  timing  error  only  at  node  D.  When  the 
resistance  of  the  shorted  path  is  5  ohms,  the  ouqtut  is  inverted  (Figure  4.8(b)).  These  faults  are  all  delected  by 
both  SPICE  and  FACT.  RSIM,  however,  predicted  an  X  value  at  node  D  when  the  fault  was  present  Because  X 
is  interpreted  as  an  indeterminate  state  rather  than  an  intennediate  voltage,  there  is  no  guarantee  that  the  X  value 
will  produce  an  effect  different  than  the  state  of  the  node  in  the  good  circuit. 

Figure  4.9  shows  a  dynamic  AND  decoder  with  a  resistive  open.  The  effect  of  the  fault  is  seen  in  Figure 
4.10  as  a  delayed  output  value.  Figure  4.11  shows  the  CMOS  latch  with  various  line-open  faults.  The  waveforms 
produced  by  SPICE  in  Figure  4.12  show  that  faults  12,  f5,  f6,  and  f7  are  detected  by  the  applied  input  patterns. 
Although  not  able  to  produce  the  detailed  waveform  such  as  SPICE,  FACT  is  accurate  enough  to  detect  all  the 
faults  detected  by  SPICE  as  seen  in  Figure  4.13.  This  shows  that  FACT  predicts  the  behavior  of  the  circuit  under 
faults  accurately  in  terms  of  voltage  levels  and  timing. 

These  simulation  results  clearly  indicate  that  the  proposed  timing  model  gives  a  significant  accuracy  to  both 
fault-free  and  fault  simulations  at  a  reasonable  cost. 
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Figure  4.8.  (a)  Waveforms  at  node  D  when  R  =  200  ohms 


Waveforms  produced  by  FACT  time  (nsec) 


Figure  4.8.  (b)  Waveforms  at  node  D  when  R  =  5  ohms 
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CHAPTER  5. 

CONCLUSIONS 


In  this  thesis,  we  presented  a  timing  model  for  fault  simulation.  The  significant  features  of  the  model  are  as 
follows. 

1)  It  uses  five  logic  values  for  more  accurate  fault-free  and  fault  simulations. 

2)  The  delay  calculation  can  handle  multiple  source  networks  and  is  therefore  a  significant  improvement  over 
the  conventional  single  source  models  of  Penfield-Rubinstein  [17]. 

3)  The  fault  model  is  very  realistic  and  incorporates  resistive  shorts  and  opens. 

4)  It  can  handle  multiple  faults  and  will  be  used  in  the  future  in  a  concurrent  fault  simulator. 

The  simulation  results  obtained  by  FACT  indicate  that  the  model  gives  more  accurate  prediction  than  that  of 
RSIM,  and  has  a  clear  advantage  over  fault  simulators  implemented  at  the  switch  level  such  as  FMOSSIM  in 
terms  of  accuracy,  all  without  much  added  cost. 
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